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I. INTRODUCTION 

The masses and properties of the SU{3) scalar mesons are long standing puzzles in 
hadron-nuclear physics related to the underlying chiral symmetry of QCD. It is also very 
interesting to describe the properties of these SU{3) scalar and pseudoscalar mesons at finite 
temperature and density. To this end it is important to study the Ul{3) x t//j(3) symmetric 
linear sigma model in the non-perturbative Gaussian Functional Approximation [l|, l2|. The 
linear sigma model is a strongly interacting renormalizable quantum field theory; due to 
the size of the self-interaction coupling constant (s) the perturbative approximations seem 
to be inapplicable. Therefore, a non-perturbative approximation, such as the Gaussian 
functional one, that is equivalent to the resummation of certain infinite classes of Feynman 



diagrams 



, l2|, is called for. 

The chiral Ul{3) x f//j(3) symmetry in the SU (3) linear sigma model is both spontaneously 
and explicitly broken, which means that some pseudoscalar mesons are Nambu-Goldstone 
(NG) bosons, i.e. with vanishing masses in the chiral limit. Straightforward solutions to the 
gap equations in the Gaussian wave functional approximation yield non-zero meson masses 
even in the chiral limit hj, ^1, however, the proof of the NG theorem used to be an open 

nti 

problem for over 30 years pL], |3|. The first solution to this problem in the 0(2) symmetric 
.gn,a .nodel was based on the Bethe-Salpete. equation 14| , but othe. p.oofs soon followed fl. 

The first proof was straightforwardly extended to 0(4) ~ SUl{2) x SUji{2) in Ref. [6| and 
was finally proven in the general 0{N) case in Ref. [3|. The SU{3) linear sigma model 
corresponds to a subgroup of the broken 0(18) symmetry, the specifics of which depend on 
the (symmetry breaking) parameters of the model, and thus readily fit into this framework, 
but the Nambu-Goldstone theorem has never been explicitly verified in the various limits of 
the chiral Ul{3) x Ur{3) Lagrangian. 

There are also several influential studies of the thermal properties of various spinless 
mesons that are based on the Gaussian approximation, both in the two-fiavor SU{2) 
and the three-flavor SU{3) cases [IQJ, but again without taking into account the Bethe- 
Salpeter equation. Therefore these studies do not obey the NG theorem in the chiral limit 
and as such are ill-suited for the study of chiral symmetry restoration. 

As for the SU{3) case, there are many studies of the properties of U{3) mesons in the 



mean field, or the Born approximation UJ, and in the Gaussian approximation 10 |. 



3 



The extension of the chiral SU{2) model to the chiral SU{3) model is not trivial, because 
there are several different self-interaction terms (three rather than one in the simplest SU{2) 
case, but one of them, the A2, is generally expected to be (much) smaller than the first one 
Ai Il3|), so we have to develop the necessary mathematical tools to deal with the scatterinsr 

nnn 

Bethe-Salpeter equation (T-matrix) for the U{3) mesons [111. Il2l. Il4j|. 

It is important to explicitly work out the NG bosons for various cases of the chiral 
Lagrangian, so as to verify which pseudoscalar mesons are NG bosons, before applying 
this formalism to non-zero temperature and/or density. The question of Ua{1) symmetry 
breaking also looms large over this endeavor, so we pay special attention to the flavor singlet- 
octet mixing. 

As this is a complicated method applied to a difficult problem, and many missteps have 
been made in the past, we take a step-by-step approach. We look first at the chiral limit: 
even here there are some non-trivial cases, such as when A2 = and c=0, (many) new naively 
unexpected Nambu-Goldstone bosons appear beyond the "elementary fields" that already 
exist in the Lagrangian - they are "composite" (bound state) NG bosons that correspond to 
the broken 0(18) symmetry rather than the Ul{3) x Ur{3) one that has (at most) nine NG 
bosons. This formation of composite NG bosons demonstrates the non-perturbative nature 
and the respect of the underlying symmetries by the GFA method. We then turn on explicit 
chiral symmetry breaking term 7^ 0, but with good SU{3) symmetry. We show how this 
explicit symmetry breaking term infiuences the masses of the pseudoscalar mesons to lowest 
(linear) approximation. 

In this paper, we present the necessary mathematical expressions necessary for the ap- 
plication of the Gaussian functional approximation, defined in Sect. Illlt to the SU{3) linear 
sigma model introduced in Sect. [TTl In Sect. [IV]we provide the expressions for the Bethe- 
Salpeter equations in various channels using the SU{3) projection operators developed in 
Sect. |Vl In Sects. IVII we verify explicitly the NG theorem for various cases and identify which 
are the NG bosons. In Sect. IVIII we briefiy discuss the role of explicit symmetry breaking 
terms by taking the simplest case. Sect. IVIIII is devoted to a summary of this paper. 
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II. THE SU{3) LINEAR SIGMA MODEL 

To understand the masses of scalar and pseudoscalar mesons, we employ the SU{3) 
linear sigma model and use the Gaussian Functional Approximation (GFA). In this section, 
we briefly review the SU{3) linear sigma model and work out the mass gap equations in the 
mean field approximation. 

The Lagrangian density of the Ul{3) x Ur{3) linear sigma model is given by 

£($) = Tr{d^,^^d''^ - m^*"^^) - Ai[Tr($t$)]2 - A2Tr($t$)^ 

+c[Det(<l>) + Det($t)] + Tr[iJ(<l> + $+)] . (1) 

The meson field matrix $ is a complex 3x3 matrix of the scalar and pseudo-scalar meson 
nonets, 

<^ = Ta(j)a = Ta{aa + tTTa) , (2) 

where cxa are the scalar fields and tTq are the pseudoscalar fields. = Aa/2 are the generators 
of U{3), where are the Gell-Mann matrices with Aq = -y/fl- The 3x3 matrix H breaks 
the chiral symmetry explicitly and is chosen as 

H = Taha , (3) 

where ha are nine (external) SU{3) symmetry breaking parameters. Only three (diagonal) 
ones, a = (0, 3, 8), are relevant and the two, a = (0, 8), are the dominant ones. In this paper, 
we only study the case ho ^ 0, and so S'f/(3) symmetry is conserved. We need to know at 



least the order of magnitude of the coupling constants. Here we may use Ref. [lOj results as 
a (rough) guide to the expected values of the coupling constants: to first approximation we 
expect Ai ~ 50, A2 — 1.5, and if we define c = Aa/^r, we find A3 ~ 50. Thus we see that this 
is indeed a strongly coupled system and that we need a non-perturbative approximation. 
The generators of f/(3) satisfy the (anti) commutation relations: 

[Aa, Aft] = 2ifabcK , and {Aa, Xb} = 2dabcK ■ (4) 
where dabc and fabc "structure constants" are defined to contain the index. The values of 



SU{3) structure constants are provided in any good textbook and in review articles 16 |. 



Those / structure constants with a zero among its a, 6, c indices are zero and those d structure 



constants containing in its a, 6, c indices are only non-zero for doab — i/S/S^ab with a,b — 
1,...,8. 

By inserting the $ field into the Lagrangian, the following Lagrangian is obtained. 

-^J^abcd{<^aO'b<^cO'd + a^^ b^^ c^^ d) + Ko'a- (5) 

The coefficients Gabc, J^abcd and Habcd are given by 

c 3 9 

Gabc = -^[dabc ~ ^i^O'^dobc + ^bodaOc + ^codabo) + 2*^000^00^60^00] , (6) 
•^ofecd = ^(^a6<^cd + ^ad^bc + ^ac^bd) + -^{dabndncd + dadndnbc + dacndnbd) , (7) 
'Habcd -^^ab^cd + -^{dabn dncd ~l~ facnfnbd ~l~ fbcnfnad^ ■ (8) 

Considering the shift of the vacuum expectation values of aa — d'a + cr'g_, the Lagrangian can 
be written as 

4 _ 4 _ 

+ (^abc - ■^J^abcd(^d)<7a<yb<^c " 3(^afcc + ■^'Habcd(^d)'^a'n'b<yc ^ '^'Habcdf^aC^bT^cT^d 
--^^abcd{<^a<7b<7c<^d + TTaTTftTTcTrd) - [/(ct) , (9) 

where we have just written Ua instead of for simplicity of writing. The potential term U {a) 
is the tree-approximation potential and aa is determined at the tree level. The tree-level 
potential is 

Tfl^ 1 
U{^a) = ~^^a ~ Gabc^a^b^c + -^^abcd^a^b^d^d " ^a^a ■ (10) 

The mean field aa is obtained by the variation. 

dU{aa) 2- o/- -- ^^x- h n (^^\ 

— — = m (Ja - 6yabc(7b(7c + -'Tabcd'^b'^c'^d " /^a = U , (iij 

5(Ta 3 

and the "tree level" masses of the scalar and pseudoscalar mesons are given by 

('mDab = rn^Sab - QGabc^c + ^^abcd^cO^d , 

imp)ab = m^Sab + QGabc^c + '^Uabcd^cO'd ■ (12) 

In general SU (3) symmetry breaking case, the 0-8 off-diagonal components of these matri- 
ces are not zero. The physical states must have a diagonal mass matrix, and we have to 
diagonalize the mass matrices. 
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III. GAUSSIAN FUNCTIONAL APPROXIMATION 

The Gaussian functional approximation (GFA) is the method based on assumine; 

the ground state solution is a Gaussian functional around the mean field. We get the effective 
potential by acting the Hamiltonian on the ground state with scalar mesons having vacuum 
expectation values. 

First the Schrodinger equation in the functional formalism is given as 

H\0) = E\0) , (13) 

where H is the total Hamiltonian and E is the corresponding energy for the wave function 
|0). The effective potential can be obtained by 

E= [ d^x {0\n\0) , (14) 



where the Hamiltonian density is obtained through the Legendre transformation as 
niaa,na) = -^^ + ^(Va,)2-iA^ + i(V7r,)2 + im2(a^ + 7r2) 

-Qabc{(^a(^b<^c - 37r„7r;,0-c) + 2'HabcdO' 

ibcd {(Ja(Jbac(Jd + TraTTbTTcTTd) - kaCTa ■ (15) 

and the ground state wave functional as a Gaussian function is 

|0) = iVexp[-i(aa - aa)G-l{m^){ab - ab) - ^7^aG-l{m^)nb] . (16) 
Here, is the normalization factor. The mass propagator is written as 

G.,.,„4../|l-^e*-.. (17) 



Finally the effective potential can be calculated as 



= (0|^|0) (18) 
= lrn'^l + ^{G-,\m,) + G-l{m^)} 

+ ^("^^ - mljGab{m^) + ^{m^ - mljGab{m^) 

- Qabci^a^b^c + "i^aiGbcim^) " Gbc{m^))} 

+ '2'Habcd{Gabifn^)^c^d + Gab{m^)Gcdim^)} 

+ -^J^abcdi^a^b^c^d + G^a^bGcdim^) 

+ 3Gab{m^)Gcd{m^) + 3Gab{m^)Gcd{m^)} 

- ha&a. (19) 
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The gap equations for scalar and pseudoscalar mesons are obtained by applying the vari- 

drua 



ational principle with respect to meson masses, -tt^ = and then the masses are given 



as 



{ml)ab = m^5ab - QGabcO-c + abcd<y c<y d + ^J^abcdG cdijn^) + ^HabcdG cdijn^) , 

("^p)a6 = m^^ab + ^Qabc^c + '^'Habcd^c^d + "i^-abcdG cdifU^) + 4:J'abcdG cdifU^) ■ (20) 

The equation following from the variation of the energy density with respect to the mean 
field value ^ = is given by 

ha = rn^aa - ^QabA^b^c + Gbc{m^) - Gbc{m^)] 

+ ^T-Labcd^bGcdim^) + ^J^abcdl'i^b^c^d + l2abGcd{rn„)] ■ (21) 

These equations provide the masses and the mean field values of the meson fields. 

IV. T-MATRIX FOR MESON-MESON SCATTERING 



As explained in Refs. 



one should work out the T-matrix for the determination of 
the pseudoscalar mesons in order to fulfill the NG theorem. We work out the a-ir scattering 
for pseudoscalar mesons. The interaction kernel in the a-n channel is written as 

4 i 4 

- iVacbd = -'i'2'Hacbd ■ 2 ■ 2 - i3{Qbea + -'Hbeaf^f)'^ ^(-03(^e(ic + -'Hedcg&g)'2 (22) 

4 _ 1 4 _ 

= —ilST-Lacbd + 36(^;,ea + -'Hbeaf^'f) :^{Qedc + -'Hedcg'^g)] ■ 

5 s — S 
With this interaction kernel we can get the T-matrix as 

-iTabcd = -iVabcd-iVabefi'n.efi-iVefcd) + ■■■ (23) 
= — '^(Kfecd + Kfee/He/Te/cd + ...) . 

Therefore, what we need to solve is the scattering matrix 

Tabcd = Vabcd + Vabef^efTefcd ■ (24) 

The polarization term for meson masses rria and rrib is 

[k - pY - ml + IE - ml + IE {2Txp 
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To work out the polarization function Hab, first let us work out the = case. In this case, 
we can write 



1 



f^2 ^2 _|_ fc2 _ ^2 _|_ ^5 y jji2 _ j^2 (^2t[)^ 



d'k loiml) - JoK) 



Here, we can write the integral as 



1 



mi — mf 



d^k 



k^ — k"^ — m? + ie 



d^k 



(26) 



(27) 



{ko - VWT^ + le) {ko + Vk^ + m? + le) (2^)' 



d'^k 



47r2 



k'^dk 



87r2 



X3 W 1 + X3 - log |X3 + \ 1 + x\ 



where X3 = A/m. We take A ~ 1 GeV to be fixed as a parameter of the model. We may 
take the 4 dimensional cut-off by transforming the above integral as fco = ik^. 



d^k 



k^ -m^ + ie (27r)4 
k^dk 1 



d'^k 



F + m2 (27r)4 



(28) 



-m 



2 [x^-logll + x^l] , 



27r2 7 k"^ + m? Att^ 
where = A/m. 

We write here the case for rria = rrib = m, which is written in the paper of Nakamura et 
al. jsj. The general case has to be worked our in the same frame. 



naa(s) = n,,(o) + 



{Any 



'1 + Jaa{s)) 



where 



JnJs) 



Am? I s 
1 arcsm 



4m2 ' 



(29) 



(30) 



for 4^ < 1, and 



Jaa{s) 



1 - 



Am^ 



log 



+ 



1 



(31) 



Am^ V 

for 1 < < 00. We should work out the general case ma 7^ mt, in the "dispersive" form 
(see Eqs. ([33])) 

r , , ■ . d'k 1 



(27r)4 - M2 + ie] [{k - P f - ^j? + le] 
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/mm(0) - ——-:KM^,{s 



1 



1 



/i 



(47r 

2 



2Ao \iJ.'^-My (47r) 



2Ao VAi^ - 



IGtt'^ / t — s — ie 



Im KMf,{t) 



where s = and the real and imaginary parts are 
Im Kmi,{s) 



-Im / 

s 



Mm I 



Re KM^,{s 



^ _ (M-^)n / _ (M + /i) 



[.-(M + /i)2) 



1 _ (^1^^ - ^^^Vanh-^ -j^ + t 

s ) Vs-{M-^i)\ 



(32) 



(33) 



V. 5C/(3) PROJECTION OPERATORS AND MIXING OPERATORS 

To solve the scattering equation (1241) . we use the method of projection operators. The 
SU (3) group structure for the SU (3) sigma model is: 

{If © Sf) ® (If © 8^) = l(i) © 8,(i) © 8j,(i) © (34) 

1(8) © 8s(8) © 8a(8) © 27s(8) © 10a(8) © T0a{8) , 



We can write out the corresponding projection operators with some manipulations: 



abed 



SnbS, 



P 



Sx{l) 
abed 



abed 



ab'Jedl ; 
a,b,e,d=0 



3 

n=l 

2 ^ ^ (^^abnd, 



n=l 



a.c=0,6,d=l---8 



a,e=l---S,,b,d=0 



P. 



1(8) 



p: 



abed 



,85(8) 
abed 



-SabS, 



abOedl 



I a,b,e,d=l---S 



^cdn 



abed 



3,275(8) 
abed 



n=l 

1 ^ 

2 ^ ^ (yfabnfedn^ i 
n=l 



a,6,c,(i=l---8 



P 



1(8) _ p85(8) 
abed abed ' 



(35) 
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p(10+10)A{8) _ - / r e e r 

^abcd - o \ "cicObd - OadObc 



_ p8A(8) 

a,b,c,d=l---a 



1 

2 

where the indices a, b, c, d can be ■ ■ ■ 8. When they are 0, they are related to li? of Eq. fl34|) : 
while when they are 1 ■ ■ - 8, they are related to Sp of Eq. (IMI) . These projection operators 
satisfy 

PabefPefcd — ^abcd ■ (36) 

We have used various relations for the derivation of the projection operators: 

(37) 



faij fbij 


— 35ab , 


daij d\jij 


— gfab , 


daij fbij 


= 0, 


daij ^ij 


= 0, 


daijdbjkdcki 


"^dabc 1 


faij fbjkdcki 


3 

2 obc ) 


faij fbjk fcki 


_ 2. 

2 J abc 1 


faij dbj k dcki 


= 0, 


dabndcdn dg^QYidbdn ~\~ dadndbcn 


= -{Sabred + SaeSbd + SadSbc) 


f acnfbdn ~l~ f adn f ben '^dg^f^^df.fjfi 


= —{^ac^bd + ^ad^bc) + ^ab^ed 


facndbdn f adn f ben fabnfedn 


= 0, 


daendbdn dg^^ndben fabnfedn 


= — ('^'^'^'^^^ ~ ^ae^bd) ■ 



Besides these projection operators, we also have several mixing operators. They are used 
to express the mixing between the singlet and octet mesons. We note that they are not 
projection operators, and so we use O to denote them. In the singlet channel, there are two 
operators: 

(38) 



^l(Ml) _ . r 
^abed - (^abOcd 



a,b=l-%,e,d=0 



^1(M2) _ . . 
'^abed — <JabOcd 



a,b=0,e,d=l-8 

They provide the mixing of the singlet and octet mesons in the resultant singlet channel. 
While in the octet channel, there are six operators 

8 

^abed — dabndcdn , , „ ' (3^) 

a=0,b,c,d=l---^ 



n=l 
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o 



8(M2) 
abed 



o 



8(M3) 
abed 



o 



8(M4) 
abed 



o 



8(M5) 
abed 



o 



8(M6) 
abed 



^ ^ dabnd, 
n=l 

8 

^ ^ dabnd, 
n=l 
8 

^ ^ dabnd, 
n=l 
8 

^ ^ dabnd, 
n=l 
8 

^ ^ dabnd. 



edn 



edn 



edn 



edn 



n=l 



a,b,d=l---8,e=0 



b=0,a,e,d=l---8 



a,b,e=l---8,d=0 



a,d=0,b,e=l- 



b,c=0,a,d=h--8 



They provide the mixing of the singlet-octet mesons and the octet mesons in the resultant 
octet channel. 

We modify the scattering equation to be 

Tabed = Vabed + Vabef^efTefed (40) 
= Vabed + Vaba'b'^a'b'e'd'Tc'd'ed , 

where Habed = ^ac^bd^ab- Then by using the projection operators as well as the mixing 
operators, we find that there is a unique expansion for Vabed of the scattering of avr — t- crvr: 



V 



abed 



^m^abed + ^S^W^abed + ^SvW^abed + 



(41) 



^^l(8)-^a6cd + ^SS{8)^abcd + '^275(8) -f^focd + ^8A{8)^abcd + *^10A(8) -f^fccd 
+ VlCAfl) f^afccd + ^nM2)<^abcd + '^8(M1) <^„6cd + ^^8(Af2) t^abcd 



,85(8) 



275(8) 



,8A(8) 



,(10+10)yl(8) 



+ Vs(M3)0 



8(A/3) 



8(M3)^abed 



+ Vs 



9'{Mi)0abcd 



8{MA) 



+ V, 



8(A-/5)C'„^^^ 



abed 
8{Mb) 



+ v. 



mm^abed 



8(A/6) 



and so does Ilabcd- Moreover, we find that for Ilabcd only the projection operators are enough, 
which means: 

ni(Mi) = ^8(Mi) = . (42) 

Since we find that both Vabed and Habed can be expanded by the projection operators fl35|) 
and the mixing operators (!38|) and (!39|) , we assume that the scattering matrix Tabed can also 
be expanded by these operators: 



rri _ rri pl(l) , rri p8x(l) p8?/(l) , 

J abed - -l-m^abed + -'■Sxil) ^abed + ^Ml)^abed + 



(43) 
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T pl(8) I T pS'^(^) -L T p27S(8) ^ p8A(8) ^ p(10+10)A(8) 

^ 1(8) ^abcd + ^ 85(8) -^afecd + ^ 275(8) -T^bcd + ^ 8A(8) -^^bcd + ^ 10A(8) -^afocd 

+ J l(Afl) f^abcd + ^l(A/2)<^afecd + ^ 8(A/1) <^afecd + ^ 8(A/2) t^'abcd 

I ^ ^8(M3) ^ ^8(Af4) ^ ^8(Af5) ^ ^8(A/6) 
+ J8(M3)<^abcd + ^8( A/4) f^afecd + ^ 8(A/5) <^afecd + ^ 8(A/6) t^'abcd ' 



and then we can separate the scattering equation fHOj) into several equations in different 
channels. We discuss them in the following subsections. 



A. Singlet Channel 



We have the following relations for all the operators: 



— n 

^abef^efcd — ^ 5 



(44) 



where denotes P^^^\ P^^^\ (9i(A^«) and 0^^*^^\ and denotes operators of other flavors. 
Therefore, we can write out the equation fHUl) in the singlet channels: 



l(M2)^aba'b' I + 



^ m^abcd + ^ m-^abcd + ^ l(Ml)<^a;,cd + ^ 1{M2)^ abed 
(^'/l(l)^afecd + ^m^bcd + '/l(A/l)<^afea'6' ' 

l^'^l(l)^afea'b' + '^l(8)-^a6a'fe' + ^HMl)<^aba'b' + '^UA/S) <^a6a'6' 

TT d1(1) TT d1(8) a 

^^l(l)-^a'b'c'd' + '-'-m^a'b'c'd' I ^ 



X 



i^^l(l)-^c'd'cd + ^l(8)-^c'd'cd + ^l(A^l)^c'd'cd + '''Hl^m'-^c'd'cd J 



It can be simplified to: 

Ti(i) 2v^Ti(jv/2) 
2\/2Ti(Mi) Ti(8) 

Vi(i) 2a/2\/i(m2) 
2\/2Vi(Mi) Vi(s) 

And its solution is 



+ 




^i(i) 2a/2Vi(m2) 

2v^\4(A/1) Vi(8) 




Ti(i) 2V2Ti(m2) 
2v^Ti(Mi) T'^s) 



Ti(i) 2a/2Ti(m2) 
2a/2Ti(mi) Ti(8) 

Vi(i) 2a/2Fi(m2) 
2V2V1(mi) ^1(8) 




(45) 



^i(i) 2A/2Vl(Af2) 

2v^\4(A'/l) Vi(8) 
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B. Octet Channel (8^.(i), 8j^(i) and 85(8)) 



We have the following relations for all the operators: 



^abef^efcd 



(46) 



where denotes p8x(i)^ p8y{i) ^ p8S{8) ^^^^ qs{Mi) ^ ^^^^ qx (denotes other operators. There- 
fore, we can similarly write out the scattering equation ( HOl) in the octet channels. After 
some simplifications it turns to be: 



8x{l) 



2rp 



— I8{M1) 



\ 



|^8(M6) 8j/(l) 
— ^-1-8(1^41) -3--'8(M4) 



Vsxil) 
|^8{M6) 



— J8(Af3) 
^85(8) J 

— VsiMl) ^ 



Vgxil) 
|V8(M6) 



i^8(M5) 



10 T/ 

— V8{A'I1) 



\ 



v 



VIOt/ 



IOt/ 

— V8{M4) 



VIOt/ 



85(8) J 



\ 



IOt/ 

— V/8(M1) 



|^8(M5) 
^83/(1) ^V8(M3) 
^85(8) 



IOt/ 

— v^8(M4) 



/ 



n 



8x{l) 






n 



-83/(1) 

n 






85(8) ) 



Tsx(l) 
3-'8(M6) 



10m 

— J8(M1) 



3^8(M5) 
^8s,(l) 
-^■l'8{Mi) 



— J8(M1) 



VlOm 
— J8(M3) 

^85(8) 



And its solution is 
f T 

3^8(Af6) 



3J8(M5) 



— ^8(M1) ^ 



T> 



83/(1) 



-3--'8(Ml) -3--'8(M4) 



VTOrp 



3 -8(M3) 
"^85(8) 



V8a;(l) 



^8(M5) -^V8(A/1) ^ 



|V8(M6) 



^8y(l) 
VTOt/ v^T/ 

-^y8^M^) -^y8^ 



— ^ 8(Af4) 



'^8(A'/1) 

2fr \/T0 



10 T/ 



fl^8(Af5) "^^8(A/1) 



|V8(A^6) ^8j/(l) -^V'8(M3) 

^ "^^8(Afl) -^V'8(M4) ^85(8) y 



3 •'8(M3) 
8S(8) 
\ 



-1 



n 



8a:(l) 






n 



83/(1) 

n 






■85(8) ) 



X 



C. 8a(8), (10 10)^(8) and 27^(8) Channels 



We have the following relations for all the operators: 



(47) 
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where denotes other operators. Therefore, we can write out the scattering equation ( HOj) 
in the 8^(8) channeh 

T8A(8) = VsA{8) + V8A(8)n8A(8)^8A(8) • (48) 

We have the following relations for all the operators: 

/'rr"™OJ/o. = 0. (49) 

where denotes operators of other flavors. Therefore, we can write out the scattering 
equation (HOi) in the decuplet channel: 



Tl0A{8) = ^10A(8) + ^10A(8)ni0A(8)7l0A(8) ■ (50) 

We have the following relations for all the operators: 

PlTO!j^, = 0, (51) 

where denotes operators of other flavors. Therefore, we can write out the scattering 
equation ( 140|) in the 27^(8) channel: 

^27S(8) = ^^275(8) + ^^275(8) 11275(8)7275(8) . (52) 

VI. THE NAMBU-GOLDSTONE THEOREM 

To check whether there are Nambu-Goldstone bosons, we need to check whether the 
scattering matrix T{s) has a pole at s = 0. In this section, we study the Nambu-Goldstone 
theorem, and verify the Nambu-Goldstone bosons in the SU{3) linear sigma model. We 
assume that there is a SU{3) symmetry, which means that only (Tq is non-zero ((Xj = 0, for 
i = 1 ■ • - 8). To further simplify our calculation, we only study the pseudoscalar channels 
(cTTT — )■ (ttt), where the pseudoscalar mesons propagate in the interaction kernel. 

Since the calculations in this system is still not so easy, our analysis will be done step by 
step. First we assume c = A2 = 0, and Ai 7^ 0. In this case, we find that all the pseudoscalar 
mesons are Nambu-Goldstone bosons. Then we assume A2 7^ 0, and find that only one 
singlet and one octet pseudoscalar mesons remain Nambu-Goldstone bosons. Finally we 
assume c 7^ 0, which is the most general case conserving SU{3) symmetry. We find that 
only one octet pseudoscalar mesons remain to be Nambu-Goldstone bosons. 
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A. Case I: c = A2 = 

This is the first step. When c = A2 = 0, we find that there is no mixing between two 
fiavor singlet mesons and among the four fiavor octet mesons in the pseudoscalar channel, 
and we can expand the potential matrix V (ct-tt scattering) by only using the non-mixing 
projection operators (1351) : 

Va,cd = V^i(i)Pf i + Vs^wPl^P + Vsya)PZi^ + (53) 

T/ pl(8) _L T/ p8-^(8) I T/ p275(8) . . pSA{8) . . p(10+T0)A(8) 

where the coefficients Vi are calculated to be 

^i(i) = 2Ai + ^fl^! , (54) 
s — [rripjoo 

^^^^^^ - + . - (mlh ' 

^82/(1) = ^1(8) = ^85(8) = ^8A(8) = ^^275(8) = ^10A(8) = 2Ai , 

To expand Ilafe at the point s = 0, first we write out its four- index form following Eq. ( I26ll : 

UaUs = 0) = 6ac6MUabis = 0) = 5,,5,,^^f^-%-ii2 , (55) 

[^s)a-a — [rripjbb 

And Habcd can also be expanded by using the projection operators (!35|) : 

n.,.. = ni(i)Pf ] + Us.^,)P'j}p + n8,(,)PS) + (56) 

n pl(S) -L n P'^'^(S) -U n p27S{8) ^ p8A(8) ^ p(10+T0)A(8) 

^^m-^abcd + ^^85(8)-r;bcd + ^h7S{8)^abcd + ^^8A(8)-f^fecd + ^^10A(8)-f^f,cd 

and its solution is: 

n8a;(l) 

n8j,(i) 

iil(8) - 1185(8) - ll8A(8) - 11275(8) - lllOA(8) - 7—3^ 7ZJ\~ ' 

In the case of c = and A2 = 0, the mass equations (120!) are diagonal, so we have 

(m|)oo = + SAi^o^ + 3XMmf) + 8Ai/o(m*J) + Ai/o(mO?) + 8Ai/o(m^^) , (58) 
{ml)u = + Xial + XJoimf) + 10Ai/o(m*^) + Ai/o(m°?) + 8Ai/o(m*^) , 









- {m],)oo 


/o(moo) 


- Io{mf) 


("^l)oo 


- {ml,)ii 


/o(m^) 


- /o«0) 


{ml)ii ■ 


- ("^p)oo 
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and 

(m^)oo = + Ai^o^ + XJoi^s') + 8Ai/o(mj^) + 3Ai/o(m°?) + 8Ai/o(m*;) , (59) 
{ml)u = + Xial + Ai/o(m°°) + 8Ai/o(mj^) + Ai/o(m°?) + IOAi/qK;) , 

where i = 1, . . . ,8, and denote (m^^)^, (m|^)^, etc.. Here we note that, as a matter of 

fact at the tree approximation level, the octet scalar mesons, the singlet pseudoscalar scalar 
meson and the octet pseudoscalar mesons all have the same mass: 

{ml)ii = (mp)oo = ("^p)n = + Ai^q • (60) 

We note that one possible consistent solution to the gap Eqs. (1211) is that they all (still) 
have the same mass: 

{ml)u = (m|)oo = {ml)u = + X^a^ + Ai/o(mO°) + 19Ai/o(mj^) . (61) 

This result is very interesting: whereas we expected to find pseudo-scalar Nambu-Goldstone 
bosons, which have zero masses, we found that even the pseudo-scalar mesons here have 
non-zero masses (mp)oo = (mp)jj. This is the usual "problem" of the NG theorem in the 
Gaussian approximation. 
From Eq. fl^Tl) . we have 

= -Xial - 3Ai/o(m°°) - 8Ai/o(m*^) - Ai/o(m°?) - 8Ai/o(mj^) . (62) 

Then using this equation together with Eqs. (158|1 and (159|1 . we obtain 

(m|)oo = 2X,al , (63) 

{ml)u = -2Ai/o(m°°) + 2Ai/o(m*j), 

{ml)oo = -2Ai/o(m°°) + 2Ai/o(m«'), 

{ml)u = -2Ai/o(m°°) + 2Ai/o(mj^). 

In order to check whether there are Nambu-Goldstone bosons, we only need to check if 
the following equations hold 

F,(s = 0)n,(s = 0) = l, (64) 

because in this case there is no flavor mixing of T-matrix elements/scattering operators. If 
Eq. holds, then the T-matrix elements subject to the following Bethe-Salpeter equation 

Ti{s) = V,{s) + V,{s)U,{s)T,{s) , (65) 
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have a pole at s = 0. This means there are massless mesons propagating, and thus the 
Nambu-Goldstone bosons turn up. 

Eq. ( I64l) can be easily checked when c = A2 = 0. By using Eqs. (1631) and (1571) . we have 



n8x(i) 



(m 



pjii 



2 ^.. ' 



and 



^8y(l) = ni{8) = n85(8) = n8A{8) = 1^275(8) = nioA{8) = , (67) 



together with Eqs. (I5lj) we have 

Vi{s = 0)n,(s = 0) = 1 , (68) 

for all the allowed flavor representations l(i), 8a;(i), 8^(1), 1(8), 85(8), 8a{8), (10©10)a(8), 275'(8). 

As there are several flavor singlets and octets, it is not clear just how many NG bosons 
in these channels are independent? Yet, it is clear that there are at least 1+8+10+27=46 
distinct NG bosons when c = and A2 = 0. That is (much) more than 9 NG bosons expected 
in the general f/L(3) x Ur{3) linear sigma model, and more than 17 NG bosons when the 
0(18) symmetry is broken down to 0(17). The explanation for the fact that there are more 
than 17 NG bosons is that the 0(18) may be dynamically broken down to a symmetry that 
is lower than 0(17), e.g. the 0(16) or even 0(15). In this sense the GFA approximation is 
substantially different from the Born, or the one-loop approximations, which are not known 
to lead to ground state(s) with "exotically broken" symmetry. 

Thus, we have proved that all the expected pseudoscalar mesons are Nambu-Goldstone 
bosons when c = and A2 = 0, but also that there are many more. As there are no 
"elementary" meson flelds in the flavor (10 © 10)^(8) and 275(8)-plets in the SU{3) linear 
sigma model, we must conclude that these NG bosons are (zero mass) bound states of 
(massive) elementary boson flelds. This goes to show that the GFA method is well and truly 
non-perturbative and capable of dynamically producing bound states even in exotic flavor 
channels, such as the (10 © 10)a(8) and 27s(8)- Of course, this does not mean that in the 
ground state of QCD there are exotic NG bosons, because the c = and A2 = conditions 
do not correspond to reality. Therefore we discuss the c = and A2 7^ case next. 
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B. Case II: c = and A2 7^ 



When c = and A2 7^ 0, the mixing between two singlet pseudoscalar mesons and 
among three octet pseudoscalar mesons exist, and we need to use the mixing operators. The 
scattering matrix V can be expanded by using the non- mixing projection operators ( I35ll as 
well as the mixing operators (!38|) and (!39|) . and the solution is 



1(1 



8x{l 
^1(8 



2 , 4 (3Ai + A2)%^ 

3 9s- {mf,)oo 



(69) 



9 s 
4 



pjii 



2, 4 A^a^ 

2Ai + -A2 + 7-4 

3 9 s — (m 



2Ai 



2, 32 

-A2 H 

3 9 s 



("^p)oo 





^85(8) 


3 9s- (mpjjj 




^8A(8) 


= 2Ai + 6A2 , 




^275(8) 


= 2A1 + 2A2, 




^10A{8) 


= 2Ai, 


^l(Ml) 


= ^l(A//2) 


2^ , 4 (3A1A2 + A2)a2 
3 9s- {mf,)oo 


V8(M1) 


= V8(M2) 


_ , 1 2 (3A1A2 + A2)a2 
3 s - (m|,)ii 


^8(M3) 


= V8(M4) 


_ . 1 2 A^ao^ 

3s- {mp)ii ' 


V8(M5) 


= V8(Af6) 


_ ^ 1 2 (3A1A2 + A2)a2 
3 s - (m|,)ii 



We do the same procedure for Hated, and the results are (after choosing s = 0) 



n 



1(1) 



("^|)c 



(70) 



n 



8a;(l) 



Io{mf) - /o(m^) 



("^l)oo 



8?/(l) 



/o(m^J) - Io{m' 
{mDii - (m|>) 



00 



n 



1(8) 



n, 



85(8) 



n, 



8A(8) 



27S{8) 



n 



10A{8) 



/o(mg) - Io{m%) 



(m 



S)ii 



[m 



pjii 



n 



l{Mi) 



8(Mi) 
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So Habcd is still "diagonal". 

In the case of c = and A2 7^ 0, the masses are diagonal, and from Eqs. fl20|) we have 



("^l)oo 



+ SXia'^Q + As^^ + (3Ai + A2)/o(m°°) + (8A1 + 8A2)/o(m; 



+ (Ai + -A2)/o(m^") + (8A1 + -A2)/o(m' 



pj 5 



(71) 



and 



S)ii 



{ml) 



00 



m' + Aia^ + A2ao + (Ai + A2)/o(m^") + (lOAi + 5A2)/oKs) 
+(Ai + ^A2)/o(mO?) + (8A1 + ^A2)/o , 



+ AiCTo' + l\2al + (Ai + ^A2)/oKO) + (8A1 + ^A2)/oK 
,00 



+(3Ai + A2)/o W") + (8A1 + 8A2)/oKi^) , 

11 17 
{ml)ii = + Xial + -Aaa^ + (Ai + -A2)/o(m°°) + (8A1 + —\2)Ioim'^] 

+ (Ai + A2)/o(m°?) + (lOAi + 5A2)/o(m^) • 



From Eq. ( 12T|) . we have 



(72) 



m 



Aido' - ^A2ao' - (3Ai + A2)/o(m°°) - (8A1 + 8A2)/o(m*^ 



(73) 



-(Ai + -A2)/o(m 



00 \ 
p ) 



(8A1 + -A2)/o(m 



Then using this equation together with Eqs. (I7T1) and (1721) . we obtain 



("^l)oo 



2Ai(Jo + TT-^S^o 



-2 



(74) 



9 1 fi 

{mlU = 2Ai(/o(m??) - /o(m°°)) + -X^m'^) - Io{mf)) + -XM^'p) - , 

{mlh = 2Ai(/o(m*;) - Io{ml')) + ^A2(/o(m«') - /o(m°°)) + ^A2(/o(m^) - /o(m^)) . 



Singlet Channel 



Since the mixing exists, we need to use Eq. 0461) derived earlier. After inserting the 
expressions of the masses and the polarization energies, Vi and Ilj, which are listed in 
Eqs. fl7^ . fl69l) and fITOl) . we can verify that at the kinematical point s = 0, we have 



2a/2Fi(mi) 



272^1(^2) \ I ni(i) \ 

^1(8) / V ° ^1(8) ; 







(75) 
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with the meaning that the determinant of the matrix within the vertical bars is zero. There- 
fore, there are Nambu-Goldstone bosons. Since there is a mixing between the two singlet 
pseudoscalar mesons, Ti(8), Tk^mi) and Ti(m2) all have a pole at s = 0. However, we 

can verify that only one of the two eigenvalues is 0. This means that only one of the singlet 
pseudoscalar meson is a Nambu-Goldstone boson. 



2. Octet Channel 8^(1) and 8s(8)) 



We calculate the solution Eq. ( H71) when c = and A2 7^ 0, and find that at the point 
s = 0, we have 



1 - 



Vsxil) |V8(M5) 



/lO 

3 



VsiMl) 



\ 



fV8(M6) 



y ^^8(M1) ^ V'8(Af4) 1^85(8) J 



8|/{1) 



— Vs(M3) 
V,: 



( 



n, 



V 



8x'(l) 










n, 



8y(i) 

n 








8S(8) j 



(76) 



where, again the vertical bars denote the determinant of the (3x3) matrix within. Therefore, 
these equations describe massless Nambu-Goldstone bosons. Moreover, we can verify that 
only one of the three eigenvalues is zero. This is difficult to prove analytically even with the 
aid of algebraic manipulation programs. Therefore, we randomly choose the values for the 
relevant parameters (coupling constants), and confirm this result. So we obtain the result 
that only one octet and one singlet of pseudoscalar mesons are Nambu-Goldstone bosons 
in this case. That agrees with the conventional result in the Born approximation, although 



the pseudoscalar spectral functions in the GFA 
simple Dirac delta function, see e.g. jisl . 



171] contain (much) more structure than a 



3. 8^(8), (10 10)^(8) and 27 s{8) Channels 

For the 8^(8) channel, we have 

+ (/o(m°°)) - /o(m«')) - 10(/o(m-j) - Ioim%)) 

For the decuplet channel, we have 

al + (/o(mOO)) - /o(m«')) - (/o(mtj) - h{m%)) 



1 — ^10^(8) nioA(8) — ^ 



(3Ai - A2)(/o(m^) - I,{m%)) + + I,{m^^) - Io{m%^)) 
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For the 27s(8) channel, we have 



1 — Vo7Q(a)Il 



27 S (8)'- '-273(8) 



A, 



cr2 + (loimf)) - JoKpO)) - 4(/o(mg) - JoKp)) 



' (3Ai - A2)(/o(m^^) - Io{mf)) + X^icx^ + Io{mf) - Io{m%^)) 
Therefore, none of them are Nambu-Goldstone bosons, as expected. 



C. Case III: c / / Aa 

In this case we assume that both c and A2 are nonzero. This is the most general case 
that conserves the SUl{3) x SU{3)fi symmetry. The scattering matrix T can be expanded 
by using the non-mixing projection operators ( 135|1 as well as the mixing operators (138|) and 
f p9|) . and its solution is 

(77) 





^l(l) 


= 2Ai 


1 2^ 1 2 2(3Ai + A2)'cr2 + 2y6(3Ai + A2)cao + 3c2 
3^9 s- (m2,)oo 






= 2Ai 


2 1 8(3Ai + X2)^ai - 4v^(3Ai + A2)cao + 3c^ 

3 ' 18 s-{ml)ii 






= 2Ai 


2 ^ 1 SA^ag - 4\/6A2Cd-o + 3c^ 

3 ^ 18 s - {mj,)ii 




Vm 


= 2Ai 


2 ^ 4 SX^a^ - 4y6A2cao + 
3^9 s - (m|,)oo 




^85(8) 


= 2Ai 


4^ ^ 5 2A^ag + 2V6A2cao + 3c^ 
3^9 s - {rr?p)ii 




^8A(8) 


= 2Ai 


+ 6A2, 




^275(8) 


= 2Ai 


+ 2A2, 




^10A(8) 


= 2Ai 




Vi(Ml) 


= Vi(A,/2) 


= ?A, 

3 


1 1 4(3AiA2 + Xl^al - y6(3Ai - A2)cffo - 
9 s- (mp)oo 


V8(M1) 


= V8(M2) 


= AsH 


1 4(3AiA2 + XI)gI + y6(6Ai + A2)ccro - 30^ 
6 s - (m\)ii 



^8(M3) — ^(M4 
V8(M5) = V8(M6 



1 4A^a2 + VeAscao - 3c2 
-^2 + 7 — 9^ ; 



A2 + 



1 8(3AiA2 + Xl)a^ - 2v/6(3Ai + 2A2)cao + 3c^ 
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s — [m 



p)ii 



We go through the same procedure for liahi and the results are (after setting s = 0) 

/o(mOO) - /o(mO?) 



n 



1(1) 



("^|)oo - ijn?p)oo 



(78) 
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n8x(i) 




- Io{m%) 


( 111 ^ ^ r^r^ 

\I"'S) 00 




^8y{l) 


Jnfmo") - 


- /o(mp ) 


{ml)ii - 


- ("^p)oo 


ni(8) 


= ^83(8) = 


n8A(8) = 


^l('mix) 


^8x(mix) 


^^8y{rnix 



/o(m*^) -/o(m^ 



which is the same as the previous case. 

In the case of c 7^ and A2 7^ 0, the masses are still diagonal, and from Eqs. fl2U]) we have 

(m|)oo = + 3Aia2 + X^a^ + (3Ai + A2)/o(mO°) + (8A1 + 8A2)/o(m*^) (79) 

+ (Ai + ^A2)/o(m°?) + (8A1 + ^A2)/o(m*;) - ^cao , 
{ml)u = + X,al + X^a^, + (Ai + A2)/oKO) + (lOAi + 5A2)/o(m|) 

+ (Ai + ^A2)/oKpO) + (8A1 + ^A2)/o(mj^) + W^cao , 



and 



(m|)oo = + Xicxl + ^2^71 + (Ai + ^A2)/oK°) + (§^1 + ^A2)/o(mf ) (80) 

+ (3Ai + A2)/o(m??) + (8A1 + 8A2)/o(m^) + ^cctq , 
{ml)u = + X^al + ^A2ao' + (Ai + ^X2)Ioimf) + (8A1 + ^A2)/o(m1^) 
+(Ai + A2)/o(m^°) + (lOAi + 5A2)/o(m^) - y^cao . 
From Eq. f l2T]) . we have 

= -Xicrl - ^X2al + ^ca^ - (3Ai + A2)/o(m^°) - (8A1 + 8A2)/o(m^^) (81) 
-(Ai + iA2)/o(mO?) - (8A1 + ^A2)/o(m^;) 

+ ^(/oK°) - Io{m'^)) - 2^-Ulo{m%) ' Io{m%)) , 
Therefore, by using this equation together with Eqs. (179|) and (IHOl) . we obtain 

(m|)oo = 2X,al + 1x2^1 - ^l\ca, + ^/f |(/o(m°°) - /o(m??)) - 2y||(Jo(mj^) - Um%)) 
{ml)u = 2Ai(/o(m^J) - /o(m°0)) + ^ca^ + '^X2al + 3A2(/o(m^^) - /o(m|)) 
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(m 



p)ii 



+ .A-(/oK°) - /o(m^°)) - 2 J|-(/oKJ) - JoK;)) , 

V 6 (To V 3 (To 

2Ai(/o(m??) - /oK°)) + y|cao + ^A2(/o(m°?) - /oK°)) + yA2(/o(m^) 

+ J\-{h{m'i) - /o(m«')) - 2./|^(/o(m^) - JoK;)) , 

V 6 (To V 3 (To 

2Ai(/o(mj^) - /o(m°°)) + ^A2(/o(m??) - I,{mf)) + ^A2(/oK;) - U^'s)) 
+ J\-{h{mf) - /o(m^°)) - 2./|^(/oKJ) - JoK;)) • 

V 0(To V 3 (To 



We calculate the solution Eq. 
s = 0, we have 



Singlet Channel 



when c 7^ and A2 7^ 0, and find that at the point 



Vi{i) 2\/2Vi(M2) 
2a/2Vi(mi) Vi{8) 




7^0 



(82) 



where, again the vertical bars denote the determinant of the (2x2) matrix within. Therefore, 
we have verified that the singlet pseudoscalar meson is not a Nambu-Goldstone boson any 
more, due to the Ua{^) symmetry breaking interaction constant c 7^ 0. 

2. Octet Channel fS^(i), ^y{i) o,nd8s(s)) 



We calculate the solution Eq. ( 147|) when c 7^ and A2 7^ 0, and find that at the point 
s = 0, we have 



Vsx(l) 



f V8(M5) -^V'8(Afl) ^ ^ 
^8y{l) 

^85(8) y 

where, again the vertical bars denote the determinant of the (3x3) matrix within. This is 
again difficult to prove analytically, so we again randomly choose the values for the relevant 
parameters, and obtain this result. Moreover, we can verify that only one of its three 
eigenvalues is zero. So we obtain the expected, yet non-trivial result that only one octet of 
pseudoscalar mesons are Nambu-Goldstone bosons. 



|^(M6) 



IOt/ 

— Vs{M3) 



VIOt/ VIOt/ 
y —V8{M1) —VsiMA) 



n8x.(l) 
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VII. EXPLICITLY BROKEN CHIRAL SYMMETRY AND DASHEN'S 

FORMULA 



When the chiral symmetry is exphcitly broken, the NG theorem turns into a relation 
between the chiral symmetry breaking parameter and the NG boson mass, as first discussed 



by D ashen 



ready been addressed in the 

n n 

% |5[. Here, we turn to the 



18lll9|. The NG theorem in the chiral limit has al 
Gaussian approximation and equivalent formalisms in Refs. 
non-chiral case. 

As shown in Ref. j4| in the chiral limit the Nambu-Goldstone particle appears as a zero- 
mass pole in the T-matrix in the pseudo-scalar channel. Next we look at the zero CM energy 
s = polarization function \4- (0)11^(0) in the non-chiral case ho = e 0. For simplicity's 
sake we only study this in the A2 = c = case (so as not to have to deal with complications 
associated with channel mixing(s) in the fiavor-singlet and octet channels). Now the gap 



Eq. (I2T1) becomes 
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We work out the BS equation for the flavor-singlet and octet channels. The polarization 
function is worked out in the fiavor-singlet channel l(i) as 
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and we obtain the similar result for the flavor-octet channel 8 



^8x(i)(o)n8,(i)(o) 



("^l)oo 



as well as the positive-parity channel flavor-octet consisting of two scalar mesons. Since 
Kr(0)n^(0) ~ 1 — we see that the pole in the s— channel propagator has moved away 
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from zero momentum. In order to find the mass, we must take into account the residue at 
the pole; thus we find = ^ + ^i^"^)^ J^st as in the Born approximation. Here Hq = e is 
the explicit symmetry breaking parameter and the "pion" symbol vr denotes the complete 
set of 17 pseudo-NG bosons, whose mass should be small if this linearized approximation is 
to hold. This result is valid only for "small" values of the explicit chiral symmetry breaking 
parameters, such as that of the SU{2)l x SUf({2) symmetry breaking that is responsible for 
the pion's mass. 

Of course, with /iq 7^ it is possible to have an explicit breaking of the 0(18) symmetry 
down to the explicitly conserved, yet spontaneously broken 0(17) symmetry. In that case 
there will remain several massless NG bosons. For instance, in the other fiavor-singlet 
channel 1(8) made up of flavor-octet mesons, we have 

\/i(8)(0)ni(8)(0) = 1, (88) 

and the same result holds for the 8j/(i), 85(8), 8^(8), 275(8) and (10 © 10)yi(8) channels. Of 
course, once one turns on A2 7^ and/or c 7^ all of these NG bosons acquire masses, as 
derived in Sect. IVI B[ The NG mesons also acquire masses when one explicitly breaks the 
0(17), SU{3) or SU{2) symmetries, e.g. by including hs ^ and/or 7^ 0. These masses 
can be evaluated by means of Dashen's formula so long as the explicit symmetry breaking 
is small, which is not the case for realistic values of /is 7^ 0; c 7^ 0. Therefore, this result is 
practically useful only for the (iso-triplet) pion masses, but not for the kaons and the r] and 
rj' mesons. This is perhaps as far as one can go using only analytic methods. 

The next step, to be taken in our next paper, will be to numerically solve the gap and 
Bethe-Salpeter equations with an explicitly broken SU (3) symmetry, so as to reproduce the 
experimental pseudo-scalar masses and their weak decay constants and thus to fix all of the 
free parameters in this model in the Gaussian approximation. 

VIII. CONCLUSION 

We have studied the Nambu-Goldstone (NG) theorem for the pseudo-scalar mesons in the 
U{3)l X Ufi{3) linear sigma model. We have constructed the ground state wave function in 
the Gaussian Functional Approximation (GFA). At this level, all the scalar mesons and the 
pseudo-scalar mesons acquire finite masses by the minimal spontaneous symmetry breaking 
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(7o 7^ 0. Hence, the NG theorem is not satisfied at the GFA level. 

Hence, we have developed a method to work out the Bethe-Salpeter equation for the 
scattering T-matrix of mesons. To this end, it was important to work out the projection 
operators in order to separate various SU{3) channels in the T-matrix. We have then 
explicitly worked out the BS equations for the pseudo-scalar mesons in the general case of 
the U{3)l X Ur{3) linear sigma model. Since the verification of the NG theorem is quite 
complicated and tricky, we have decided to work out the NG theorem step by step. 

Firstly we studied the Ai 7^ and X2 — C — O case. In this case, we verified that the NG 
bosons appear in the usual fiavor-nonet channel, where the NG bosons are present at the 
mean field, or the Born approximation level. Additionally, we have found new composite 
NG bosons in certain other flavor channels that correspond to the breaking of the extended 
0(18) symmetry down to a lower symmetry. 

Then we studied the case with Ai 7^ and A2 7^ 0, but c — 0. In this case, we found the 
usual fiavor-nonet of NG bosons. We then studied the case when all the coupling constant 
in the Lagrangian are non-zero. In this case, we found only the fiavor-octet pseudo-scalar 
mesons as the NG bosons: the ninth pseudo-scalar meson acquires a non-zero mass and thus 
is not an NG boson any more. Of course, c 7^ corresponds to the explicit U^(l) symmetry 
breaking, that affects the rj and rj' mesons, and is comparable with, or perhaps even larger 
than the exphcit breaking of the SU{3)l x SUr{3) symmetry. 

We have discussed another simple case in order to examine how low-mass pseudo-NG 
bosons emerge due to the explicit chiral symmetry breaking: when the Lagrangian has just 
one small explicit chiral symmetry breaking parameter ho = s 0. There we confirmed 
that Dashen's result for pseudo-NG boson masses hold in the Gaussian approximation. This 
result is valid only for "small" values of the explicit chiral symmetry breaking parameters, 
such as that of the SU{2)l x SUr{2) symmetry breaking that is responsible for the pion's 
mass. 

In this paper, we have analyzed the appearance of NG bosons for various cases in the 
chiral U{3)l x Uii{3) linear sigma model Lagrangian. We have also studied the effect of the 
"small" explicit chiral symmetry breaking term to provide a small mass to the pseudo-scalar 
bosons. This result is practically useful only for the (iso-triplet of) pions, but not for the 
kaons and the 77 and 77' mesons. 

This is perhaps as far as one can possibly go using only analytic methods. The next 
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step, to be taken in our next paper, will be to numerically solve the gap and Bethe-Salpeter 
equations with an explicitly broken SU{3) symmetry, so as to reproduce the experimental 
pseudo-scalar masses and their weak decay constants and thus to fix all of the free parameters 
in this model in the Gaussian approximation. Then it will be possible and (very) interesting 
to calculate the spectra of scalar bosons in the SU{3) chiral sigma model. 
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